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Abstract 

We show that partons traversing a quark-gluon plasma can lose substantial amounts 
of energy also by scatterings, and not only through medium-induced radiation as 
mainly considered previously. Results from Monte Carlo simulations of soft inter¬ 
actions of partons, emerging from a hard scattering, through multiple elastic scat¬ 
terings on gluons in an expanding relativistic plasma show a sizeable jet quenching 
which can account for a substantial part of the effect observed in RHIC data. 


Key words: jet quenching, quark-gluon plasma, QCD, scattering 
PACS: 24.85.-Fp, 12.38.Mh 


Partons emerging from hard scattering processes and traversing a quark-gluon 
plasma (QGP) are expected to lose a substantial amount of their energy 
through interactions with the plasma, resulting in a suppression of high-p^ 
particles [1,2]. Such jet quenching has indeed been observed at RHIC [3,4,5,6] 
and is believed to be not only a signal for QGP formation, but also a tool for 
investigating the properties of the plasma. So far theoretical efforts have con¬ 
centrated on medium-induced gluon bremsstrahlung as energy loss mechanism 
[7,8,9,10,11,12] based on arguments that the collisional energy loss is small 
[13,14,15]. In this letter, however, we demonstrate that energy loss caused by 
multiple elastic scatterings is not negligible, but can contribute signihcantly to 
the observed jet quenching effect. The importance was also realised in [16,17] 
and scattering was included at the same level as radiation in [18]. 

Such multiple scatterings can be related to recent developments to understand 
soft QCD interactions of a high energy parton with a colour background held. 
This has, in particular, been used as a novel way to understand diffractive hard 


Preprint submitted to Elsevier Science 


2 February 2008 




scattering in ep or pp collisions [19]. In the Soft Colour Interaction (SCI) model 
[20,21] a hard-scattered parton interacts with the proton remnant via soft 
gluon exchanges. Here, it is the exchange of colour charge which is important, 
since this changes the colour topology of the event such that hadronisation 
will produce a different hnal state, e.g. diffractive ones with a gap in the 
rapidity distribution of hadrons. The phenomenological success of the SCI 
model indicates that it captures the most essential QCD dynamics for soft 
gluon exchanges. A theoretical basis for this model has recently been found 
[22] in terms of QCD rescatterings of an energetic parton with the target 
colour held via one or more gluons as expressed through the Wilson line in 
the parton density dehnition. 

In this letter, we develop the SCI model to apply for a parton that rescatters in 
a quark-gluon plasma. The essential quantity for the energy loss through such 
elastic scattering is then the energy-momentum transfer involved, rather than 
the colour exchange. Since the QGP is much denser and of larger volume than 
the colour background held of a single proton, the parton should experience 
many more interactions so that even a small momentum transfer per scattering 
could add up and result in a sizeable energy loss. 

The SCI jet quenching model [23] is implemented as a Monte Carlo event 
generator based on Pythia (version 6.2) [24] which is used to simulate hard 
scattering processes based on perturbative QCD 2 —> 2 matrix elements and 
initial and hnal state parton showers based on DGLAP evolution. These per¬ 
turbative QCD (pQCD) processes are not altered, but are used as in standard 
high-p^ pp collisions. However, before the emerging partons hadronise, we 
introduce the possibility that they may interact with the plasma. 

This interaction of the hard parton with the plasma is treated as elastic scat¬ 
tering on the partons in the QGP, with a squared momentum transfer t. Al¬ 
though this is treated as a 2 —> 2 parton-parton scattering (in the cm frame), 
we cannot use pQCD matrix elements since that would only apply for the small 
cross-section processes with large t. Instead we want to examine energy loss 
through the dominant soft interactions where perturbation theory does not 
apply. Being unable to calculate a proper t-distribution from hrst principles, 
we assume that it can be approximated by a Gaussian distribution having a 
width at as a parameter of the model. Naturally, the mean momentum transfer 
should not exceed a few hundred MeV in order to be non-perturbative. 

The partons emerging from the hard pQCD phase are traced through the 
plasma and scatter, with a probability Pint, on each plasma gluon along their 
way within a screening radius Pscr- We use Pint — 0.5 since this was found as 
a universal parameter value for the SCI model to ht rapidity gap data from 
both ep and pp collisions [20,21], as well as charmonium formation [25,26]. 
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The QGP is represented by an ideal relativistic gluon gas, while quarks are 
neglected because the QGP is initially gluon rich and the gluons come much 
faster into thermal equilibrium than the quarks. This means that the number 
and energy densities of gluons are connected to the temperature via Ug = 
ggT^2 ({3)/2Tr‘^ and eg = vr^^f^T^/SO so that Ug oc 

The time evolution is governed by a Bjorken-like longitudinal expansion [27], 
i.e. with the equation of state of an ultra-relativistic ideal gas the time de¬ 
pendence of the energy density and temperature is given by e(r) oc and 
T(r) oc where r = is the proper time. Therefore, the density 

of gluons drops very fast, namely as ng(r) oc r“^. It is assumed that the local 
initial energy density is proportional to the number of binary nucleon-nucleon 
collisions at impact parameter b, i.e. e{x,y,b) oc T\u{x — b/2,y) ■T\u{x + b/2,y) 
where the nuclear thickness function Tau is estimated with a simple Glauber 
model [28]. 

The parameter that governs the energy density of the plasma is Cq, the energy 
density in most central collisions (i.e. 6 = 0) at Tq = 1 fm/c averaged over 
the transverse area. This fixes the normalisation of the density profile for any 
centrality. 

The centrality of a nucleus-nucleus collision is defined by the fraction of the 
total geometrical cross section it takes. A centrality class, i.e. a range in cen¬ 
trality can be translated into an impact parameter range using the Glauber 
model calculation. In the simulation an impact parameter is chosen for each 
event in a given range of centrality according to da oc bdb (ignoring the 
fluctuations in the experimental quantity used to determine the centrality). 

The hard scatterings are treated as in pp collisions and are distributed in the 
transverse plane according to the number of binary nucleon-nucleon collisions 
per unit transverse area, which is obtained from the Glauber model calcula¬ 
tion. As discussed above, the emerging hard partons from the pQGD processes 
are then tracked through the plasma. It is assumed that no interactions occur 
before the formation of the QGP at r,, but thereafter the local gluon density is 
updated in each time step taking the changes due to expansion and the energy 
density prohle into account. The procedure is stopped when either the parton 
under consideration leaves the QGP or the local temperature drops below Tc, 
the critical temperature of the phase transition. 

Hadronisation of partons emerging from the plasma presents new interesting 
problems. The conventional models developed for e’''e“ annihilation and ap¬ 
plied for ep and pp, need not be applicable but there is little guidance yet how 
the presence of a quark-gluon plasma affects the fragmentation of energetic 
partons. It is also desirable to disentangle the effects of direct energy loss and 
modihed hadronisation. Therefore, in a hrst step, the standard fragmentation 
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procedure is used in this model. However, while the standard Lund string frag¬ 
mentation model [29] is the best option to use for the pp reference, it is not 
easily applied to the case of partons emerging from a plasma after several soft 
colour exchanges which have changed the colour topology. It is quite unclear 
where the string from such a parton should be connected and if the concept 
of the normal string applies at all. The pragmatic way out that we have taken 
in this first study is to instead apply independent hadronisation [30,31] for 
our simulations of AuAu collisions, but keeping the same basic fragmenta¬ 
tion function f{z) for the energy-momentum fraction 2 ; given to the produced 
hadron in each iteration. 

Finally, one also has to account for the Cronin effect, i.e. the p± broadening of 
the hnal state hadrons due to conventional initial state scatterings. We have 
included this according to the model suggested in [32,33], i.e. the variance of 
the intrinsic /c^-distribution is increased by a constant a for each scattering 
prior to the hard interaction, i.e. crl^{x, y, h) = q -1- a ■ {Nsca.t{x, y, b) — 1). 
In our model does not depend on the hard scattering momentum transfer 
scale because parton showers are treated explicitly. The parameter a was 
hxed independently with the help of dAu data at the same energy. 

Table 1 


Parameters of the model for the quark-gluon plasma (using input from hydrody¬ 
namics and lattice calculations) and the SCI jet quenching mechanism. 


Parameter 

value 

obtained from 

QGP formation time 

n 

0.2 fm/c 

based on saturation scale [35] 

energy density at tq = Ifm/c 

eo 

5.5 GeV fm-3 

fixed from hydro [36] 

critical temperature 

Tc 

0.175 GeV 

fixed from lattice [34] 

gluon mass 

rUg 

0.2 GeV 

chosen here 

interaction probability 

Pint 

0.5 

fixed from SCI [20,21] 

screening radius 

^scr 

0.3 fm 

cf. [37] 

width of Gaussian t distr. 

(Tt 

0.5 GeV^ 

chosen here 

increase of per scattering 

a 

IGeV^ 

fitted from dAu data 


The parameters of the model are listed in Table 1. Not all of them are free, 
the critical temperature is taken from lattice calculations [34]. The energy 
density has been determined with hydrodynamic calculations [36] and can 
vary only in a small range when considered at a fixed time Tq = Ifm/c. The 
formation time is chosen in accordance with saturation scale considerations 
[35]. Using the known time evolution the energy density at any time can 
be calculated, in particular at the time when the QGP is formed. The 
soft interaction probability -Pint is hxed from the SCI model. The gluons in 
the plasma have an effective mass, which is a free parameter. This mass has 
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important dynamical consequences for the energy loss, since a larger mass 
results in smaller energy loss. In fact, it was the assumption of static scattering 
centres, i.e. inhnitely massive plasma partons, that gave vanishing energy loss 
through scattering so that only losses through the medium-induced radiation 
was treated in [2]. The width of the momentum transfer distribution at is free 
and may be regarded as the most important parameter because it regulates 
how much energy can be lost per collision. 


The results are discussed in terms of the nuclear modihcation factor and the 
two-particle azimuthal correlation. The nuclear modihcation factor dehned as 


Rab{p±,v) = 


N, 


evt 


. / (^^coii) dVPP \ ' 
dp±dT]) \ dp±dp) 


( 1 ) 


is a measure of deviations of the p^-spectra obtained in nucleus-nucleus col¬ 
lisions from the pp result scaled with the average number of binary collisions 
(^coii)- The two-particle azimuthal correlation 


D(A0) = 


1 


diV 


iVtrig d(A0) 


( 2 ) 


is the azimuthal separation from a trigger particle, normalised to the number 
of trigger particles in the data set. Here, both the trigger particle and the 
other particles are required to be in certain ranges of p±. 

The model results for the nuclear modihcation factor (above 3 GeV) are shown 
in Fig. 1 together with the PHENIX data [38] for diherent centrality classes. 
Also shown are the results for d-l-Au collisions that where used to determine 
the Cronin parameter a. The resulting Cronin ehect of a broadened intrinsic 
/i;_L-distribution is shown separately for the most central and the most periph¬ 
eral AuAu collisions, demonstrating that it is bigger in central events where 
the number of collisions of a nucleon is larger for geometrical reasons. 

The model reproduces the approximately hat shape of the RauAu data. Al¬ 
though the model is in agreement with the data in peripheral collisions, it 
reaches only ~ 50% of the suppression for the most central collisions. The 
comparison of the SCI jet quenching model with the data in Figure 2 illus¬ 
trates again that the model is in agreement with the data for peripheral and 
semiperipheral collisions but shows a weaker dependence on centrality when 
going to more central collisions. Nevertheless, even in the most central colli¬ 
sions, the model shows a strong jet quenching effect which should be taken 
relative to the increase in i?AuAu due to the Cronin effect. It should be noted 
in this context that interactions with the hadronic hnal state can also result 
in a signihcant suppression of high-p^ hadrons [39] and this effect would have 
to be added to any QGP suppression. 
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Fig. 1. Nuclear modification factor, eq. (1), for different centrality classes (0 — 10% 
represents the most central and 80 — 92% the most peripheral collisions) of gold-gold 
(and minimum bias deuterium-gold in lower right panel). PHENIX data [38] com¬ 
pared to the SCI jet quenching model (starting at p_L = 3 GeV due to the cut-off 
used on Pythia’s hard scattering matrix elements). The result from initial state 
scattering, but no interactions with the plasma, is shown by the curves ’Cronin 
effect only’. 

In order to reproduce the data in central collisions with the SCI jet quenching 
model alone one would have to roughly double the overall energy loss. There 
are several ways how this can be achieved: One could for instance increase the 
number of scatterings by increasing the screening radius from 0.3 fm to 0.5 fm, 
or one could double the momentum transfer by increasing at from 0.5 GeV^ to 
2 GeV^. While with the latter we think we are leaving the reasonable scale for 
soft interactions one could consider the former, which effectively increases the 
in-medium parton-parton cross section from 3 to 8 mb ^. It is not sufficient 
to increase the number of scatterings by a factor 2 because the energy loss 


^ In a recent study the parton-parton cross section had to be increased to 45 mb to 
describe elliptic flow data in a parton-cascade [40] 
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in a single scattering decreases with decreasing parton energy so that the 
parton loses less energy in the additional scatterings when its energy is already 
reduced. With either of these two parameter variations the SCI jet quenching 
model agrees with the data for the most central collisions, but the centrality 
dependence is not linear so that the calculation falls below the data in semi¬ 
central and, to a lesser degree, peripheral events (Fig. 2). 

The centrality dependence is sensitive to the energy and path length depen¬ 
dence of the energy loss mechanisms and can be used to discriminate between 
different types of model as investigated in [41]. In particular, models in which 
the energy loss is based on coherent processes and grows with the path length 
squared, seem to describe the linear centrality dependence better. However, 
also the energy dependence is relevant for the dependence of RauAu on central¬ 
ity. In the SCI jet quenching model, using our default parameters, the energy 
loss in a single scattering is negligible for small parton energies (< 1 GeV de¬ 
pending on the plasma temperature), then it rises steeply and flattens at high 
energies reaching roughly 400 MeV at 20 GeV in the case of light quarks. 



Fig. 2. The nuclear modification factor RauAu (averaged for with p± > 3 GeV) 
versus centrality for AuAu data [38] and for the SCI jet quenching model, with 
default parameters and with the screening radius Rgcr increased from 0.3 to 0.5 fm. 
The quenching effect should not be considered relative to i?AuAu = 1 (dotted line) 
corresponding to no nuclear effects, but relative to the upper curve including the 
Cronin effect. 

The resulting azimuthal correlation. Fig. 3, show that in peripheral collisions 
there is a clear jet-like peak at A0 = 0 and a somewhat lower and broader one 
at A(j) = TT, much like in pp interactions. In central collisions, however, there 
is a clearly suppressed away side jet, although the peak does not disappear 
as in the data [42,43]. To understand the quenching of the away side jet is a 
general problem and we are not aware of any successful model. Therefore, we 
discuss it here in more general terms. 

The amount of energy that a parton loses through the interactions with the 
QGP is determined by the number of interactions and the energy loss in a 
single scattering. The latter is dominated by the width of the t-distribution, 
but depends also on the gluon mass and the temperature. A higher gluon 
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Fig. 3. Model results for the 2-particle azimuthal correlation, eq. (2), between a 
trigger particle (having 4GeV < < 6GeV and defining zero azimuth, but 

not included in the plot) and other particles (having p± > 2GeV) in peripheral 
(80 — 92%) and central (0 — 10%) collisions. The model is in reasonable agreement 
with the data for peripheral collisions, but in central collisions the away-side peak 
disappears completely in the data [42]. 

mass or a higher temperature lead to a smaller energy loss (the energy of the 
glnons in the QGP is proportional to T and, dne to kinematics, the energy 
loss is largest when the glnon has a small energy). While the QGP expands 
the temperatnre drops according to T oc and the energy loss mechanism 
becomes more efficient at later stages, bnt the density of glnons in the plasma 
drops mnch faster, namely n oc so that the dilation exceeds the effect of 
cooling. Fnrthermore the amonnt of energy lost by a parton depends mnch 
stronger on the glnon density than on the temperatnre. Gonseqnently, the 
energy loss at early times dominates the total energy loss. 

Related to this is onr finding that the largest energy loss is obtained when 
the hard scattering takes place in the centre of the overlap region. This is 
illnstrated in Fignre 4 which shows the energy loss of light qnarks in a central 
collision for different emission points. When moving the hard emission point 
from the centre towards the snrface the energy loss gets smaller becanse even 
those qnarks that move towards the centre (0 = tt) cannot have a longer 
effective path length (i.e. glnon density integrated along the path) in the plasma 
than those coming from the centre dne to the expansion of the plasma. By 
the time a qnark coming from a distance (r) reaches the centre, the density 
has already dropped and, moreover, for part of the qnark’s path-length there 
will even be no plasma to interact with dne to the plasma’s limited life time 
of only 5.2fm/c compared to the ~ 14 fm diameter of the gold nnclens. 

Another point is that the locations of the hard interactions are, in fact, con¬ 
centrated towards the centre of the QGP where the nnmber of binary nncleon- 
nncleon collisions per nnit transverse area is highest. The path lengths of the 
two hard-scattered partons are thns typically similar which makes it nnlikely 
that one of them loses mnch more energy than the other (which is needed for 
the disappearance of one jet). Again, the limited plasma life time (in rela- 








tion to the gold nucleus size) prevents large asymmetries in the path lengths 
through the plasma (the path length difference cannot become larger than the 
plasma lifetime). 


These three points prevent large asymmetries and may explain why it is so 
difficult to get a substantial suppression (or disappearance) of the away-side 
jet. On the other hand, the inhomogeneous energy density distribution can 
help to amplify small asymmetries. Furthermore the hard interaction is at 
RHIC energies dominated hy q + g^q + g scatterings, which should also help 
because the gluon interacts more strongly than the quark and should thus lose 
more energy. These effects are included in our simulations, but are not strong 
enough to create a large enough asymmetry with associated strong quenching 
of the away side jet. 




Fig. 4. Energy loss AE of light quarks emerging from a hard scattering with 5 GeV 
energy at different distances r from the centre in central AuAu collisions (6 = 0) 
with azimuthal angle 4> (as illustrated to the right). Quarks emitted at large r with 
(j) ^ TT have a long path through the collision region, but by the time they reach 
the centre (where the QGP is densest) the density has already dropped to very low 
values due to the expansion or the QGP has even already hadronised. 

It should be noted that this discussion is not only valid for this particular 
model, but is of a more general nature since the arguments are independent of 
the details of the interaction mechanism. It should thus apply to all scenarios 
in which the energy loss depends strongly on the density of scattering centres, 
which is also the case for medium-induced gluon radiation. 

In fact, it seems that the two-particle azimuthal correlation is sensitive to 
details of the fragmentation procedure and is therefore afflicted with an addi¬ 
tional uncertainty. Already Lund string fragmentation and independent frag¬ 
mentation with the same fragmentation function lead to quite different asso¬ 
ciated multiplicities. A softer fragmentation function produces more hadrons 
but with lower momentum and since the parton p± spectrum is steeply falling 
already small changes in the fragmentation function might cause that the 
hadrons fall below the p± threshold of 2 GeV used for the azimuthal correla¬ 
tions. 
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In conclusion, energy loss due to soft scattering of energetic partons in the 
QGP can contribute significantly to the jet quenching observed at RHIC. 
The SCI jet quenching model gives a nuclear modihcation factor having the 
correct p± dependence and a magnitude which can account for most of the 
effect observed in peripheral collisions and about half the effect in central 
collisions. There is also a suppression of the away-side jet, although not as 
strong as in data. This depends on the distribution of hard scattering events 
and initial energy density as well as the plasma evolution, which are not specihc 
to our particular model. This may give handles for further investigations of the 
quark-gluon plasma. The centrality dependence of the observed jet quenching 
may indicate the need for taking into account the coherence between individual 
scatterings in the plasma. For an improved understanding of the jet quenching 
phenomenon, one needs to take into account both energy loss through medium- 
induced radiation and through scattering. 
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